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Resonance tongues are mode-locking regions of parameter space in which stable 
periodic solutions occur; they commonly occur, for example, near Neimark-Sacker bi- 
furcations. For piecewise-smooth, continuous maps these tongues typically have a dis- 
tinctive lens-chain (or sausage) shape in two-parameter bifurcation diagrams. We give 
a symbolic description of a class of "rotational" periodic solutions that display lens- 
chain structures for a general iV-dimensional map. We then unfold the codimension- 
two, shrinking point bifurcation, where the tongues have zero width. A number of 
codimension-one bifurcation curves emanate from shrinking points and we determine 
those that form tongue boundaries. 

1 Introduction 

Mathematical models often incorporate a nonsmooth component in order to describe a phys- 
ical discontinuity or sudden change. Piecewise-smooth, continuous maps are a class of such 
systems. We say that a map 



can be partitioned into countably many regions where F has a different smooth functional 
form. The map is nondifferentiable on codimension-one region boundaries called switching 
manifolds. 
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Piecewise-smooth, continuous maps have been used as models in many areas, for example, 
economics PQ[2], power electronics [21 HI [2], and cellular neural networks [5]. Furthermore, 
they arise as Poincare maps of piecewise-smooth systems of differential equations, particu- 
larly near sliding bifurcations [TJ [8] and near so-called corner collisions [9| [TO]. 

Piecewise-smooth, continuous maps may exhibit border- collision bifurcations when a fixed 
point crosses a switching manifold under smooth variation of system parameters [UJ [12], [131 
Q3]. A fundamental characteristic of these bifurcations is that multipliers associated with the 
fixed point may change discontinuously at the bifurcation. In some cases, border-collision bi- 
furcations resemble bifurcations of smooth maps (such as saddle-node and flip bifurcations); 
however, invariant sets created at the bifurcation generically exhibit a linear (not quadratic) 
growth in size under variation of system parameters. Alternatively, border-collision bifur- 
cations may have no smooth analogue (for example, an immediate transition from a fixed 
point to a period-three cycle in a one-dimensional map). Furthermore, these bifurcations 
may be extremely complicated, producing multiple attractors and chaotic dynamics. 

Resonance (or Arnold) tongues are regions in parameter space within which there is an 
attracting periodic solution. In some cases (e.g., circle maps), the periodic solutions in a 
tongue can be characterized by their rotation number. As was first observed for a one- 
dimensional, piecewise-linear, circle map [IS], two-parameter pictures of resonance tongues 
for piecewise-smooth, continuous maps typically exhibit a distinctive lens-chain (or sausage) 
geometry [T6J [T7J [TJ [18], [19], as in Fig. [TJ The boundaries of a resonance tongue correspond 
either to a loss of stability of the associated periodic solution or to a collision of one point of 
the solution with the switching manifold. The latter case usually corresponds to the collision 
of the stable periodic solution with an unstable periodic solution of the same period and is 
known as a border-collision fold bifurcation. (This need not always occur, as is shown in 
|16j.) Following [15] . we call the codimension-two points at which resonance tongues have 
zero width, shrinking points (sometimes called "waist points"). The purpose of the present 
paper is to explore the bifurcation that occurs near such points for a general N- dimensional 
map in the neighborhood of a switching manifold. 

A switching manifold locally divides phase space into two regions. Any nearby period-n 
orbit with no points on the switching manifold consists of, say, I points on one side and 
n — I points on the other side of the switching manifold. We observe that I differs by one for 
coexisting stable and unstable periodic orbits in lens shaped resonance regions. Moreover, 
the number / also differs by one for stable periodic orbits in adjoining lenses. For example, 
Fig. [2] shows a magnification of Fig. [TJ near a generic shrinking point in the 2/7-tongue. The 
borders of the lens near such a shrinking point correspond to border- collision bifurcations 
of the two orbits in the lens, and thus to the occurrence of an orbit with a point on the 
switching manifold. 

Consequently, a shrinking point is defined geometrically by the occurrence of a periodic 
solution that two points on the switching manifold. However, in §7] we show that there 
are two cases that must be treated separately: terminating and non-terminating shrinking 
points. Non-terminating shrinking points are defined algebraically, but by the singularity of 
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Figure 1: Resonance tongues corresponding to period- n cycles up to n = 30 of (T5]) with 
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6 = e\ and // = 1, as shown in [IB]. The 



rotation number associated with each resonance tongue is indicated. 



certain matrices (the 'border collision" matrices). Terminating shrinking points correspond 
to the case I = 1 or n — 1; they are also defined algebraically by the occurrence of certain 
eigenvalues. We will prove in Lem. [7] (under certain nondegeneracy conditions, see Defs. [2] 
and [3]) the following. 

• The algebraic shrinking point conditions imply the geometrical condition of having two 
points on the switching manifold. 

• At a shrinking point, the map has an invariant polygon (that is typically nonplanar) 
on which the dynamics are conjugate to a rigid rotation. 

Since a shrinking point is defined by two conditions (two points on the switching mani- 
fold), the unfolding of a shrinking point bifurcation requires the variation of two parameters. 
In the generic two-parameter picture, a (non-terminating) shrinking point appears, at first 
glance, to lie at the intersection of two smooth curves. However, this is an illusion: there are 
really four distinct curves corresponding to four different border-collision fold bifurcations. 
We will prove in Thm. the following. 

• Each of the four border-collision curves is quadratically tangent to one of the others 
and appears on only one side of the shrinking point. 
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• The extension of these curves through the shrinking point corresponds to virtual solu- 
tions. 

• At each border-collision fold bifurcation one particular point of a particular period-n 
cycle lies on the switching manifold. 

• There are curves in parameter space along which each point of this cycle lies on the 
switching manifold, see Fig. [2j however except for the four that form the resonance 
tongue boundaries, these correspond to virtual solutions. 




Figure 2: A magnification of Fig. Q] showing the 2/7-resonance tongue near a shrinking point. 
Sketches of four period-7 orbits are shown. The orbit with I = 3 (its symbol sequence is LLRRLRR, 
see $3|) is admissible in both gray regions. In the upper region this orbit is stable; in the lower 
region it is a saddle. Let {x{\ denote the points on this orbit. On the curves labeled (i) the point Xi 
lies on the switching manifold. By Th. [9l the curves (0), (1), (3) and (5) form tongue boundaries. 
A saddle period-7 orbit with I = 2 is admissible in the upper region and coincides with {x{\ on 
the boundaries (0) and (1). Similarly, a stable period-7 orbit with I = 4 is admissible in the lower 
region and coincides with {x{\ on the boundaries (3) and (5). At the shrinking point there is an 
invariant heptagon. On the curve V the linear solution system (|13p is singular. 

The remainder of this paper is organized as follows. In $2] we introduce the general map 
that will be studied throughout this paper and determine fixed points. Symbolic dynamics 
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are introduced in §3] to describe periodic solutions. The periodic orbit corresponding to 
a particular symbol sequence may be found by solving a linear system, see §HJ In §0 we 
define the new concept of a "rotational symbol sequence" and using symbolic dynamics are 
able to classify periodic solutions as either rotational or non-rotational. The motivation is 
that rotational periodic solutions correspond to lens-chain resonance regions whereas non- 
rotational periodic solutions may not [16]. In §|6]we discuss the number of distinct rotational 
periodic solutions. Formal definitions for terminating and non-terminating shrinking points 
are provided in §7J This section describes the singular nature of shrinking points and the 
construction of invariant, nonplanar, polygons. This mostly preparatory work sets up §S]that 
gives the unfolding of the shrinking point bifurcation summarized by Thm. [91 Concluding 
comments are given in §|9j 



2 A framework for local analysis 

It is well-known (see for instance [12]) that dynamical behavior of the piecewise-smooth, 
continuous map ([1]) local to a nondegenerate border-collision bifurcation of a fixed point at 
a smooth point on a switching manifold is described by a piecewise-affine map of the form 

T > _ f M _ [ s < o 

X - U[X) - 1 f*(x), s > ' [Z) 



where 

for i = L, R. Here 



f i Jx) = /ib + A i x, (3) 



s = ejx, (4) 

denotes the first component of x G R N , \i e R is a system parameter, b e M. N , and A\_ and 
Ar are N x N matrices, which, by continuity of (j2J), must be identical in all but possibly their 
first columns. Nonlinear terms only affect local dynamics in degenerate cases (for example 
when Ai has a multiplier on the unit circle). Also, for simplicity, we have assumed the 
switching manifold is analytic and therefore may be transformed to the plane, s = 0. 

The map ([2]) is a homeomorphism if and only if det(A\_) det(^R) > 0. Furthermore, ([21) 
has the scaling symmetry 

fx^(Xx) = Xf,(x), VAX), (5) 

and consequently every bounded invariant set of (|2J) collapses to the origin as /i — > 0. Fur- 
thermore, it suffices to consider only \i = —1, 0, 1. 

Potential fixed points of (|2J) are obtained by solving x*^ 1 ' = fl(x*^ 1 '). Whenever Ai does 
not have a multiplier of one, x*^> is unique and given by 



x 



*(i) 



yL{I-Ai)~ x b. (6) 
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Each x*^\ as given by (jSJ), is a fixed point of (jSJ) and said to be admissible whenever s*W 
has the appropriate sign (negative for % = L and positive for i = R), otherwise it is called 
virtual. We have 

where adj(X) denotes the adjugate matrix of a matrix X, (adj(X)X = det(X)J). Since A\_ 
and Ar differ in only possibly the first column, their adjugates share the same first row, 
which we denote by g T : 

g T = e]adj (/ — A\_) = e]adj (I - A R ) . (7) 

Thus 

S *« = ^ h (8) 
det(J - At) ■ 1 } 

The condition g T b ^ is a nondegeneracy condition for the border-collision bifurcation, 
guaranteeing that the distance between each x*^ and the switching manifold varies linearly 
with the bifurcation parameter /x. Furthermore, if (I — AA is nonsingular and a, denotes 
the number of real multipliers of Ai that are greater than one, then det(J — AA is positive 
if and only if is even. Thus, as described in [20], by (jHJ), if cll + clr is even, a;** 1 ) and £*( R ) 
are admissible for different signs of /x, thus as /x is varied through zero a single fixed point 
persists. Conversely if ai + a^ is odd, a;*( L ) and x*^ coexist for the same sign of /x. At /x = 
the fixed points collide and annihilate in a nonsmooth fold bifurcation. 



3 Symbolic dynamics 

We refer to any sequence, S, that has elements taken from the alphabet, {L, R}, as a symbol 
sequence. This paper focuses on periodic solutions; thus, we will assume S is finite and let 
n denote the length of S. We index the elements of S from to n — 1. Arithmetic on the 
indices of S is usually modulo n. For clarity, throughout this paper we omit "mod n" where 
it is clear modulo arithmetic is being used. 

Here we introduce notation relating to symbol sequences. Given n e N, the collection 
of all symbol sequences of length n is {L, R} n = E^. The i th left cyclic permutation is an 
operator, crj : Ti 1 ^ ~ > ^25 defined by (<JiS)j = Si + j and we use the notation, £W = <JiS. The 
i th flip permutation is an operator, \i '■ ^2 — * ^21 that flips the 2 th element of S (i.e. L —>■ R 
and R — > L) and leaves all other elements unchanged. We use the notation, S" 1 = Xi<~>- F° r 
example if S = LRLRR then S 3 = LRLLR and = LRRLR. In general &® = (S 1 )^ ^ 
(S^) 1 = S^ 1 because with the same example S 3 ^ = LLRLR and S (2 ^ = LRRRR. 

The i th multiplication permutation is an operator, 7Tj : — > ^2 > defined by (niS)j = <Sy . 
Notice (ninjS)k = {-njS) ik = S ijk = (vr ii >S) fc thus 

TtiTtj = TCij . (9) 
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Consequently 7Tj is an invertible operator if and only if gcd(i,n) = 1 and the inverse of 
7Tj is 7Tj-i (where i~ l is the multiplicative inverse of i modulo n). Also notice (<JiiTjS)k = 
(TTjS)i +k = Sij+jk = (a ij S) j k = (Ttj(JijS) h . Hence 

a-iKj = HjOij . (10) 

We let SS denote the concatenation of the symbol sequences S and S and let S k G 
Yi\ n denote the symbol sequence formed by the concatenation of k copies of S. A symbol 
sequence is called primitive if it cannot be written as a power, S k , for any k > 1. S is 
primitive if and only if S ^ for all i ^ 0. 



4 Describing and locating periodic solutions 

Each orbit of ([2]) can be coded by a symbol sequence that gives its itinerary relative to the 
switching manifold. However, instead of defining symbol sequences for orbits, we find it 
preferable to do the reverse. Given a point x = x G M N , we denote Xj as the i th iterate of 
x under the maps and fj} in the order determined by 5 G SJ: 

x i+1 = f*( Xi ) (11) 

In general this is different from iterating x under the map ([2]). However, if the sequence {x{\ 
satisfies the admissibility condition: 

s = f L, whenever Sj < . . 

1 R, whenever > 

for every i, then {xj} coincides with the forward orbit of x under (T5]). When ( TL21) holds for 
every i, {xi} is admissible, otherwise it is virtual. 

For a given symbol sequence S, we are interested in finding x G ~R N such that x = x n , 
because then {x ,Xi, . . . ,x n _i} is a period-n orbit. We call this orbit an S-cycle. iS-cycles 
are determined by the linear system 

Xi = A So x + fib , 
x 2 = A Sl xi + fib , 



xo = A Sn ^x n -i + fib . 
Elimination of the points xi, . . . , x n _i, gives 

(/ - M s )x = fiP s b . (13) 

where 

M s = A 5n _ 1 . . . Asq , (14) 
P s = I + As n _ 1 +A Sn _ 1 A Sn _ 2 + --- + A Sn _ 1 ...As 1 . (15) 
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We call (fTBl) the n-cycle solution system of S. If (J — M$) is nonsingular, then (TlBl) has the 
unique solution 

x = u(J - M s )- l P s b . (16) 

Stability of the period-n orbit is determined by Ms and for this reason we call Ms the 
stability matrix of S. In view of Lem. H] (see below), we call Ps the border- collision matrix 
of iS. Notice Ps is independent of <So, thus 

Ps = P s o . (17) 

Also, it is easily verified that Ms and M s o differ in only their first column. Consequently, 
the map that describes the n th iterate of xq under either S or <S°: 

/ fxPsb + Msx, s<0 
71 \ fiP s b + M s oX } s>0 ' 1 J 

is piecewise-smooth continuous and has the same form as (T5]). For this reason (T5]) may be used 
to investigate dynamical behavior local to border-collision bifurcations of periodic solutions. 

We now state five fundamental lemmas relating to n-cycle solution systems that we will 
utilize in £j7]and £0 Lems. H] and [5] are generalizations of those given in [16] . 

Lemma 1. Suppose x solves the n-cycle solution system (T73)) of S and Sj = ejxi = 0. Then 
x also solves the n-cycle solution system of S l . 

Proof. By continuity: A\_Xi = A^Xi, hence there is no restriction on the i th element of S. 

□ 

Lemma 2. Suppose x and x solve the n-cycle solution systems of S and S° respectively. 
Then det(J - M s )s = det(J - M s o)s. 

Proof. By flT3]) and flHD, we have (J - M s )x = {I - M s o)x. Since (J - M s ) and (I - 
M 5 o) are identical except in the first column, the first row of their adjugates are identical. 
Multiplication on the left by this row to the previous equation yields the desired result. □ 

Lemma 3. For any i, det(7 — M s(i )) = det(J — M s ). 

Proof. Suppose w.l.o.g., S = L. If A L is nonsingular, then (/ — M s m) = A L (I — M S )A'[ 1 
which verifies the result for i = 1. Since nonsingular matrices are dense in the set of all 
matrices and the determinant of a matrix is a continuous function of its elements, the result 
for i = 1 is also true even when A\_ is singular. Repetition of this argument completes the 
result for any i. □ 

Lemma 4. Suppose (I — Ms) is nonsingular, /i^O and g T b ^ 0. Then the point xq, given 
by [Tb]). lies on the switching manifold if and only if Ps is singular. 
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det(/ - M s ) ^ 


det(/-M 5 ) = 


det(P 5 ) 7^ 


unique solution 
and s ^ 


no solution 


det(P 5 ) = 


unique solution 
and So = 


possibly uncountably 
many solutions 



Table 1: A grid summarizing the nature of solutions to ([TBI when fi ^ and g T b ^ 
as determined by Lems. H] and Ms is the stability matrix of S, fTHl) . and Ps is the 
border-collision matrix of S, ffTo^) . 

Lemma 5. Suppose Ps is nonsingular, |i / and g T b ^ 0. Then the n-cycle solution 
system ( T73j) /ias a solution if and only if (I — Ms) is nonsingular. 

Proofs for Lems. H] and [5] are given in appendix [A] An interpretation of these last two 
lemmas is presented in Fig. [TJ The situation det(7 — M s ) = det(Ps) = is generically 
co dimension-two. When appropriate nondegeneracy conditions are satisfied it is equivalent 
to a shrinking point, see §0 

5 Rotational symbol sequences 

Of particular interest is the situation that the map ([2]) exhibits an invariant, topological 
circle that crosses the switching manifold at two points. When the restriction of the map 
to this circle is homeomorphic to a monotone increasing circle map we find that periodic 
solutions on the circle can only have certain symbol sequences. We call these sequences 
rotational symbol sequences. 

Definition 1 (Rotational Symbol Sequence). Let /, m, n e N, with l,m < n and gcd(m, n) = 
1. Let S = S[l,m,n} be the symbol sequence of length n defined by 

where d is the multiplicative inverse of m modulo n, i.e. dm = 1 mod n. Then we say S and 
any cyclic permutation of S is a rotational symbol sequence. 

Notice d always exists and is unique because m/n is an irreducible fraction [21]. For 
example if (/, m, n) = (3, 2, 7) then d = 4, hence So = <S>4 = <5>8 mod 7=1 = L, thus <S[3, 2, 7] = 
LLRRLRR. 

A pictorial method for computing S in terms of /, m and n is to select n points on a 
circle, then draw a vertical line through the circle such that / points lie to the left of the line, 
see Fig. |3j Label the first point to the left of the lower intersection of the circle and line, 
point 0. Move m points clockwise from and label this point 1. Continue stepping clockwise 
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labeling every to point with a number that is one greater than the previous number until 
all points are labeled. Then <Sj = L if the point % lies to the left of the vertical line and R 
otherwise. 




Figure 3: Illustration of the pictorial method for determining <S[3, 2, 7] = LLRRLRR. 

The circle represents an invariant circle of (J2J) and the vertical line represents the switching 
manifold. Taking to steps clockwise corresponds to evaluating the map (J2J) once. Thus 
S[l, m, n] is the symbol sequence of an n-cycle of ([2]) that has / points to the left of the 
switching manifold and with rotation number m/n. 

The following lemma states some basic properties of rotational symbol sequences. 

Lemma 6. Suppose S[l,m, n] is a rotational symbol sequence. 

a) S[l, m, n] = n m S[l, 1, n] . 

b) If S[l, m, ti]q = L and S[l, to, n]_ d = R then i = 0. 

c) S[l,m,n] is primitive. 

d) S[l, m, n] , n — TO, 77. 1 . 

e) 5[Z,TO,n]^- x ) d ) 5 = S[l,m,nf( ld l 

f) If < toi,TO2 < f are distinct integers coprime to n and I ^ l,n — 1, t/ien <S[i,TOi,n] 

noi a cyclic permutation of S[l,rri2,n\. 

Proof. Recall the notation o~iS = for the i th left cyclic permutation and Xi<S = <S l for 
the flip of the i element. 
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a) Since dm = 1, 

{iT m S[l, 1, n}) id = S[l, 1, n\i = L if and only if % = 0, . . . , I - 1 , 
matching the definition of «S[Z,m, n], (TE?j) . 

b) Let j = im (equivalently, z = j<i). Then S[l, m, n}^ = S[l, m, n]jd = L implies < j < 
I — 1. Similarly, S[l, m, n]_ d = S[l, m, n]y_i)d = R implies / < j — 1 < n — 1. The only 
value of j that satisfies both inequalities is zero, hence i = 0. 

c) From part (b), S[l, m, n] differs from any nontrivial cyclic permutation of itself in either 
the th or the (— d) th element. Therefore, S[l,m,n] is primitive. 

d) By definition, 

S[l, m, n]id = L if and only if i = 0, . . . , I — 1 , (20) 
and S[l, n — m, n]^ n -d) = L if and only if i — 0, ... ,1 — 1 , (21) 

since the multiplicative inverse of (n — m) is (n — d). From (1201) , by letting j — i — l + 1 
we obtain 

S[l, m, n]f~ l)d) = L if and only if j = -I + 1, . . . , , (22) 
thus S[l, m, = L if and only if j — 0, . . . , I — 1 , 

which matches ( 12 II) . 

e) From (1221) we obtain 

5[Z, m, nf d l ~ 1)m = L if and only if z = -I + 1, . . . , -1 . (23) 

Also 

S[l, m, n]° id = L if and only if i = 1, . . . ,1 — 1 , 
thus S[l, m, n]f d ld) = L if and only if j = -I + 1, . . . , -1 , (24) 

where we have set j — i — I. Equation (|24p matches (J23J) proving the result. 

f) Let d\ and c?2 denote the multiplicative inverses of mi and m 2 modulo n, respectively. 
Let 

S = ir dl S[l,mi,n] , (25) 
<S = 7T dl (S[l,m 2 ,n]M) , (26) 

where k G Z. Since 7^ is an invertible operator, it remains to show that for 
any fc 6 Z. 
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Using (jHJ) and part (a) we find 

S = S[l,l,n\. (27) 

Re-expressing S in the form we desire is a little more complicated but requires no more 
that the basic known properties concerning multiplicative permutations, tt. 

S = (7r dl S[l,m 2 ,n]Y km ^ , by P} 

= {n dun n m2tn S[l,l,n]) ikmi) , by part (a) 

= (7r Ainl S[/, 1, n]) (fcmi) , where m = d x m 2 , by © 

= <S[Z,m,n] (femi) , by part (a). (28) 

Notice m 7^ 1 since mi 7^ m,2- Also m 7^ n — 1 since, otherwise, d\m2 = — 1 =>■ 
(n — di)m 2 = 1 =>• r«2 = n — mi (since the multiplicative inverse of (n — d{) is 
(n — mi)), thus mi +m2 = n which is a contradiction since mi, m.2 < f by assumption. 

Using f|2^|) . by the definition of an rotational symbol sequence ( TT^l) . if 5j = L, then 

^j+d = L for all but one value of i & {0, . . . , n — 1} (where (im = 1). We now show this 

property of S is not exhibited by S and hence S and S cannot be equal. 

Using ( 1271) and remembering d, I ^ 1, n — 1, if d < I, then <S t _ d = S t _^ +1 = L and 

Si = Si + i = R. Similarly if d > I, we have 5o = <S>i = L and <S d = <S d+1 = R. In either 

case 

S^S (29) 

for any k G Z. 

□ 



6 The cardinality of symbol sequences 

Let Af n [A/^ ot ] denote the number of primitive symbol sequences [primitive rotational symbol 
sequences] of length n that are distinct up to cyclic permutation. To use combinatorics 
terminology, M n is the number of mbead necklaces of two colors with primitive period n, 
and it is the number of binary Lyndon words of length n [22| 123] . The formulas we give for 
AT n and Af* ot use the Mobius function: 

{1 if n = 1 or n is the product of an even number of primes 
— 1 if n is the product of an odd number of primes (30) 
otherwise, i.e. n is not square free 

and Euler's totient function: 

(f(n) = the number of positive integers less than n and coprime to n (31) 

where <^(1) = 1. The first few values of fi(n) and (p(n) are: 
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n 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 




1 


-1 


-1 





-1 


1 


-1 








1 


-1 





ip{n) 


1 


1 


2 


2 


4 


2 


6 


4 


6 


4 


10 


4 



Then 

a/; = -VmV (32) 

a|n 

A^ ot = 2 + — ^— </?(n) , forn > 3 (33) 

where Yl a \n denotes summation over all divisors, a, of n. The first formula, (132]) . is well- 
known (see for instance [221 [23] for a derivation). To arrive at (1331) we count the number 
of distinct (up to cyclic permutation) S[l,m, n] for an arbitrary fixed n: When I = 1 there 
is one distinct rotational symbol sequence, namely LR n_1 . Similarly when I — n — 1 there 
is only L n_1 R. There are (n — 3) values of I left to consider. For each of these, there are 
ip(n) possible values for m. By Lem. 00, the values of m that are less than | yield 
distinct symbol sequences and by Lem. IHEI the remaining values of m only produce cyclic 
permutations. Thus we have fl33l) . 

The first few values of Af n and J\f* ot are: 



n 


i 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


K 


2 


1 


2 


3 


6 


9 


18 


30 


56 


99 


186 


335 







1 


2 


3 


6 


5 


14 


12 


20 


16 


42 


20 



All primitive symbol sequences of length n < 6 are rotational symbol sequences except when 
n = 1 because we do not consider L and R to be rotational. There are four distinct primitive 
symbol sequences of length six that are not rotational. These are: LRLRRR, LLRLRR, LLRRLR 
and LLLRLR. Roughly, as n increases the number of distinct primitive symbol sequences of 
length n that are not rotational also increases. To quantify this statement, Af n (1321) . grows 
like e n , whereas A/"™* ( |33|) . grows like n 2 (since y?(n) grows linearly (tp(p) = p — 1 for any 
prime p)). Thus for large n the majority of primitive symbol sequences are not rotational. 

7 Shrinking points 

Roughly speaking, as in [15] , we call points where lens-chain shaped resonance tongues have 
zero width, shrinking points. The aim of this section is to provide a rigorous foundation 
for the unfolding of shrinking points described in §BJ We define two classes of shrinking 
points: terminating and non-terminating. This categorization provides a distinction be- 
tween shrinking points that lie at the end of a lens-chain (terminating) and those that lie 
in the middle (non-terminating) (recall Fig. [p. We assume associated symbol sequences are 
rotational (see §5$); this assumption is crucial to our analysis. Our main results are Lem. [71 
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which states that at a shrinking point there exists an invariant, nonplanar (though planar 
in special cases) polygon, and Cor. [HI which shows that shrinking points are a hub for the 
singularity of important matrices, 

Near a shrinking point, a lens shaped resonance tongue corresponding to the existence 
of an admissible, stable period-n cycle, {xj}, has two boundaries. These correspond to 
border-collision fold bifurcations of {xi} with an unstable orbit of the same period. At 
each boundary, one point on the orbit lies on the switching manifold. The first question to 
address is the following: which points lie on the switching manifold at the two resonance 
tongue boundaries? 

Suppose {x^ has an associated symbol sequence that is rotational, S[l,m,n\. We may 
picture the points, Xj, lying on a topological circle, as in Fig. [3j The switching manifold 
intersects the circle at two points. If this structure is maintained as parameters vary, then 
it seems reasonable that only points that lie adjacent to an intersection can collide with the 
switching manifold. When I ^ l,n— 1, there are four such points: xq, X-d, and xu- 

(For example, if [I, m, n] = [3, 2, 7), as in Fig. [31 then d — 4 and the four adjacent points are 
xo, X3, x\ and x$.) 

Suppose w.l.o.g. that Xq lies on the switching manifold at one resonance tongue boundary. 
In the interior of the tongue the symbol sequence of the corresponding unstable periodic 
solution then differs from S in the th element, i.e., is S°. If we assume that {x^} collides 
and annihilates with the same unstable periodic solution on the second boundary, it must 
be X(7_i)d that lies on the switching manifold there because the only index i for which S %d is 
a cyclic permutation of S°, is i = I — 1. Consequently, in view of Lem. 0] since {x^ is not 
always well-defined, we define non-terminating shrinking points by the singularity of P$ and 
P S ((i-i)d) (Def. [2] below). In §S]we will show that resonance tongue boundaries at which the 
remaining two points, and xid, lie on the switching manifold, form a second lens shaped 
resonance tongue emanating from the shrinking point. 

As we will see, the cases / = 1 and I = n — 1 correspond to terminating shrinking 
points. It suffices to consider only one of these cases, we choose I — n — 1, because they are 
interchangeable via swapping L and R. For the remainder of this paper we assume / 7^ 1 
which greatly simplifies our analysis. In particular we find it is always reasonable to assume 
that the unstable periodic solution described above always exists (though it may not be 
admissible). That is (/ — M s q) is nonsingular and therefore the n-cycle solution system, 
(TT3"j) . of S° has the unique solution 

p = Po = Kl-M s o)- l P s ob. (34) 

We denote the i th iterate of p via the symbol sequence S° by pi and let U denote its first 
component. The orbit, {pi}, will play a pivotal role in our analysis. 

Definition 2 (Non- Terminating Shrinking Point). 

Consider the map ([2]) with N > 2 and suppose that /i 7^ and g T b 7^ 0. 

Let S = S[l, m, n\ be a rotational symbol sequence with 1 < I < n — 1. 



14 



Suppose 

Ps and P S ((i-i)d) are singular. (the singularity condition) 

Let S = S° and S = S ld and assume (I — M§) and (/ — M£) are nonsingular. 

Suppose the orbit, {pi}, of (|34|) . is admissible. 

Then we say ([2]) is at a non-terminating shrinking point. 

Def. [2] essentially characterizes non-terminating shrinking points as the codimension- 
two phenomenon at which the border-collision matrices P$ and P s ai-i)d) are simultaneously 
singular. The following definition for terminating shrinking points is quite different because 
the defining characteristic of these points, in our opinion, is the co dimension-two requirement 
that one fixed point (here x*^ because we are assuming I — n — 1) is admissible and has a 
pair of associated multipliers on the unit circle with a particular rational angular frequency. 

Definition 3 (Terminating Shrinking Point). 

Consider the map (T5]) with iV > 2, suppose (/ — A\_) is nonsingular and det ^j_^ L ) < (i.e. the 
fixed point x*^ is admissible, see ([8])). 

Let S = S[l, m, n] be a rotational symbol sequence with I = n — 1 and n > 3. 

Let S = S° and suppose (J — M§) is nonsingular. 

Suppose 

e ± « are multipliers of A\_ . (the singularity condition) 
Then we say (j2j) is at a terminating shrinking point. 



P(i-Dd 



Pu 



<M)d 



Figure 4: A schematic illustrating iterates of p, ( 1341) . on or near the switching manifold at 
a non-terminating shrinking point. 

In two dimensions, terminating shrinking points are center bifurcations for rational ro- 
tation numbers, m/n, studied in [241 125] . These authors show there exists an invariant 
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polygon that has one side on the switching manifold, within which all points other than 
the fixed point belong to periodic orbits with rotation number, m/n. As we will show, in 
higher dimensions this behavior occurs on the center manifold of x*^ corresponding to the 

i . i. I 27rim . 

multipliers e « , call it h c . 

The following lemma concerns the orbit, {pi}, of (jMj) . The reader should take care to 
notice that singularity of the matrix Pg^-yd), by Lem. [U implies that pu (not P(/_i)<z) lies 
on the switching manifold because 

P$W = Ps°( ld ") = -^sCC'- 1 )^ = PsiV-W > (35) 
where the second equality is Lem. [(Jo and (JIT]) is used for the last equality 



A 




Figure 6: The invariant, nonplanar polygon, V, at a non-terminating shrinking point cor- 
responding to the rotational symbol sequence <S[2,2, 5] = LRRLR for the map §2§ when 
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Lemma 7. Suppose |D) a shrinking point. Then, 
a) t = t ld = 0; 

&j if the shrinking point is non-terminating, then t^, < and > as in 

Fig. [7J- z/ the shrinking point is terminating, then t id < for all i ^ 0,-1; 

c) {pi} has period n; 

d) the iterates of p are the vertices of an invariant, nonplanar n-gon, V , that is comprised 
of uncountably many S-cycles and the restriction of to V is homeomorphic to rigid 
rotation with rotation number, m/n. 
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Separate proofs for non-terminating and terminating shrinking points are given in Ap- 
pendix [Bl If the shrinking point is terminating, then the polygon V is planar. In general, 
when N > 2, V is nonplanar as in Fig. 6. In any case V has two vertices on the switching 
manifold. 

V is comprised of many 5-cycles, in other words the n-cycle solution system of S has 
more than one distinct solution. Therefore at a shrinking point the matrix (J — Mg) is 
singular. Furthermore, for each i, by Lem. 0(1 — M 5 (i)) is singular. However, there exist 
solutions to the n-cycle solution system of (such as p^, so to avoid a contradiction with 
Lem. [51 P 5 (i) must be singular. Consequently we have the following: 

Corollary 8. Suppose |D) is at a shrinking point. Then, 

a) (I — M s ) is singular, 

b) P S {i) is singular, for all i. 

The singularity of a matrix is a co dimension-one phenomenon. Thus by Cor. [HI we expect 
there will be curves in two-dimensional parameter space passing through shrinking points 
along which (I — Ms) and each P^w is singular. Some of these curves form resonance tongue 
boundaries, and it is these that form the focus of the following section. 



8 Unfolding shrinking points 

This section presents an unfolding of the dynamics near terminating and non-terminating 
shrinking points. We begin by assuming (T5]) is parameter dependent and perform a coordi- 
nate transformation such that, locally, in two dimensions, two tongue boundaries lie on the 
positive coordinate axes. 

Suppose (T5]) varies with a set of parameters £ G WL K (K > 2) and is C k+1 (k > 2). Then 
A\_ and Ar are C k functions of £. Suppose (j2J) is at a shrinking point when £ = 0, for some 
fixed (jl 7^ 0. Then, det(Ps(£)) and det(P 5 cci-i)d)(£)) are C k scalar functions that have the 
value zero when £ = 0. If the 2x2 Jacobian d ^ det ^ Ps ^ t ^s(< l - 1 ) d ) } j s nc ,nsingular for some i th 

and j th components of £, we may utilize the implicit function theorem to obtain a coordinate 
transformation of £, such that det(Ps) = when rj = & = and det(P s ai-i)d)) = when 
v = £j : = 0. In what follows below we omit dependence of (j2j) on the remaining [K — 2) 
components of £. 

Let x(t]j v) denote the unique solution to the n-cycle solution system, f[T3|) . of S = S° 
for small values of rj and [i. Let p = x(0, 0) (to coincide with (j34j) ). If the shrinking point 
is non-terminating, let x(r], v) denote the unique solution to the n-cycle solution system of 
S = S ld for small values of rj and v. Let and Sj denote the first components of Xi and X{, 
respectively. From the above assumptions on n and fi, by Lem. H] we have s(0, v) = and 
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hdiVi 0) = 0) ( see (05])). If ^-s(0, 0) and J^s^(0, 0) are both nonzero, we may redefine r] and 
v via a nonlinear scaling so that 

3(7;, i/) = 77(1 +0(1)), (36) 
s ld ( V ,v) = 1/(1+0(1)). (37) 

Note, we use 0(k) (and o(k)) to denote terms that are order k or larger (larger than order 
k) in all variables and parameters. 

Theorem 9. Suppose the map (Tj|) is a C k+1 (k > 2) function of parameters 77 and v, has a 
shrinking point when 77 = v — and ioe /icrae ( [375)) and p?7| ). Write 

det(J - M 5 (r/, i/)) = fcir; + A; 2 z/ + 0(2) , (38) 

and assume ki, ki ^ 0. v4s nsna/ tue assume N > 2. 

Let = {(77, z/) I 77, v > 0}. T/ien for small rj and v, the S and S -cycles are admissible 
in ^1 and collide in border- collision fold bifurcations at the boundaries, 77 = and v = 0. 

// the shrinking point is non-terminating, there exist C k functions gi, g 2 '■ K. — > M. such that 
for small 77 and v , the S and <S -cycles are admissible in ^ 2 = {(v-> u ) \ V — 9i{ u )i v — 92{v)} 
and collide in border- collision fold bifurcations at the boundaries, 77 = g%{v) and v = g 2 (if). 
Furthermore, the lowest-order, nonzero terms of gi and g2 are g'{(0), g'^ifi) < 0. 

Near non-terminating shrinking points, Th. [9] predicts the bifurcation set sketched in 
Fig. [6j The four resonance tongue boundaries that emanate from the shrinking point corre- 
spond to the collision of the points of the iS-cycle: x , x_ d , Xn-i)^ %id, with the switching 
manifold. As mentioned in §7] these are the four points adjacent to the switching manifold in 
the interior of the tongue. 

By applying the implicit function theorem to (I38I) there exists a C k function 

%) = -^ + 0(t/ 2 ), (39) 
k 2 

such that det(J — M s (r], h(r)))) = 0. But by Th. [9j S-cycles exist throughout the first 
quadrant of parameter space therefore the curve v = h(r]) must not enter this quadrant. 
Hence, since h'(0) = 

hk 2 <0, (40) 

as in Fig. [6j 

By Cor. [HI the border-collision matrix of every cyclic permutation of S is singular at a 
shrinking point. Th. [9] describes curves along which some of these matrices remain singular. 
We now describe all other such curves. If po and pu are the only points of the period-n 
cycle, {pi}, that lie on the switching manifold (as is generically the case), then for each 
i ^ 0,1 — 1,1, —1 if the shrinking point is non-terminating and for each % 7^ 0, —2, —1 if the 
shrinking point is terminating, there exists a C k function 

Qi(v) = -%^r; + 0(t; 2 ) , (41) 
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that satisfies det (Pg(id) (i], (?;(??))) = 0- m eacn case t(i+i)d and have the same sign, thus 
by (HDj) we have g.(0) < 0. Although the border-collision matrix, P S (id), is singular along 
v = qi(i]), these curves are not seen in resonance tongue diagrams because the related 
periodic solutions are virtual. 

So far we have described curves passing through the shrinking point along which P 5 (m) 
is singular for each i, except for i — —1 when the shrinking point is terminating. Here 
S (~d) = RL™- 1 , hence P^-d) = I + A L + ■ ■ ■ + A^ 1 = {I - A L y l {I - A^). The matrix A"l has 
a multiplier 1 with an algebraic multiplicity of at least two, therefore det(P 5 (-d) (rj, v)) = 0(2). 
Thus in this case, for small rj and z/, generically P^c-d) is singular only when r] = v = 0. 

In addition, for terminating shrinking points there generically exists a curve in parameter 
space passing through the shrinking point along which the matrix, A\_, has a pair of complex 
multipliers on the unit circle. As one moves along the curve, the angular frequency of 
multipliers changes. Whenever the angular frequency is rational there may be additional 
terminating shrinking points. Such curves are described in [161 12H [221 [T9] . 

The border-collision bifurcations occurring on resonance tongue boundaries are described 
by piecewise-affine maps of the form f|T8|) . Th. [9] states that these bifurcations are nonsmooth 
folds, thus we may use Feigen's results [2D], summarized in §2J to determine the relative 
stability of periodic solutions near shrinking points. Let a^S] denote the number of real 
multipliers of Mg that are greater than one in the interior of ^j. Then ai [S] + ai[S] and 
o 2 [5] + a 2 [S] are odd. Also, det(J - M s (v,v)) = 0(1), (PD, thus |ai[5] - a 2 [S]\ = 1. 
Consequently, if iS-cycles are stable in the interior of ^i, then iS-cycles are unstable and 
5-cycles are unstable in the interior of ^2- 




Figure 6: A schematic of the generic bifurcation set in a neighborhood of a non-terminating 
shrinking point. 
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Finally we turn to the proof of Th. [9j 
Proof of Th. Write 

det(I - M s (rj,v)) = h + 0(l), (42) 

det(J - M s {r], v)) = k 4 + 0(l), (43) 

where k 3 , fc 4 ^ by assumption. Let 

x(r], v) = n{I - M S { V , v))- x P s {n, v)b , (44) 

whenever v ^ h(rj), ([39]) . For small v ^ 0, if rj = then x is defined and coincides with x 
because det(Ps) = 0. Therefore Xi d (0, v) = x^(0, u), for all i. In particular, 

s«(0, i/) = t w + 0(u) , (45) 

so that by ([36]) and (j37j) 

s(0,i/) = 0, 

s w (0,i/) = u + 0(u 2 ). (46) 

Similarly for small r\ ^ 0, if v = then x is defined and det(P S ((!-i)d)) = 0, thus x solves 
the n-cycle solution system of S^ l ~^ d = <S^ d \ Therefore Xi d (r},0) = X(i + i) d (r], 0), for all i. In 
particular, 

Sid(v, 0) = t(i+i)d + 0(rj) , (47) 

so that 

s_ d (77,0) = r/ + 0(r/ 2 ), 
S(/_i)d(?7,0) = 0. 

We now compute lowest order terms in the Taylor series of si d (r], v). An application of Lem. [2] 
to the symbol sequence , using also Lem. [3] produces 

det(J - M s {r), v))si d (r), v) = det(I - Mg(r), v))s ld (r), v) . 

Using ([38D, dH, (BSD and g7]) we obtain 

MV,0) = ^^V + 0(v 2 ). (48) 

By the implicit function theorem, there exists a unique C k function g\ : R — ► R such that 
for small u, si d (gi(u), v) = and 

^i(^) = -t-^ — ^ 2 + 0(^ 3 ) . (49) 
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In a similar fashion, by applying Lem. [2] to the symbol sequence 5 (_c!) and noting S^ = S, 
we obtain 

s(0,u) = M^z/ + 0(z/ 2 ), (50) 
K4 

so that the implicit function theorem gives the function 

92(V) = -T^hv 2 + 0(rf) . (51) 



We now prove k\k 2 > 0. Let Ej C M 2 denote the intersection of the interior of the 



th 



quadrant with a sufficiently small neighborhood of the origin (i = 1, ... ,4). Suppose for a 
contradiction, k\k 2 < 0. Then, by (1551) . (I — Ms) is nonsingular throughout E 2 and hence s 
is continuous throughout E 2 . Notice s = only when 77 = 0, thus the sign of s is constant 
in E 2 . When z/ = 0, by admissibility and using (1471) we have s(r/, 0) = + 0(1) < 0, thus 
s is negative throughout E 2 . In particular, s is negative when v = ^ for small 77 < 0. By 
applying Lem. [2] to 5, we obtain 



det(J - M 5 (t7, v))s(r), i/) = det(/ - M s {r], i/))sfa, v) . (52) 
n 

k 2 



Substituting z/ = ^ and using ([361), (EHJ and flUS} gives 



(2A; 1 r7 + 0(r / 2 )) S (r / ,^) = (fc 3 + 0(»|))(f| + O^ 2 )) , 
k 2 

=> 2^3(7],^-) = k 3 + 0( V ), 

k 2 

=> hh < . (53) 

Via a similar argument with S(/_iv, we find fc 2 /C3 < 0. This inequality provides a contradiction 
with (1551) . thus kik 2 > 0. Hence s and S(j_i)d are in fact continuous and negative throughout 
Ex and E 3 , thus fci^ < 0. Via similar arguments we find kik 4 > 0. Consequently we have 

sgn(fci) = sgn(fc 2 ) = -sgn(fc 3 ) = sgn(/c 4 ) (54) 

Furthermore, from P9l) and (jBTj) we have 

^(0),^'(0) <0, (55) 

since t {l+1)d ,t_ d > 0. 

Near (77, z/) = (0, 0), since {p^ is admissible, 5-cycles are admissible if and only if s, su > 
0, thus only in ^ by (1561) and (157j) . Similarly 5-cycles are admissible only if s, sm < 0. By 
(150|) . §§(0, 0) = is positive by (154)) and admissibility. Thus £(77, z/) < when v < h 2 {rf). 

Similarly, by (1451) . si^(t],u) < when 77 < h\(v). Therefore 5-cycles are admissible in \l/ 2 . 
The result for iS-cycles follows by looking at the signs of each Sj on the boundaries of ^1 and 

*2- □ 
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9 Discussion 



Dynamics local to a nondegenerate border-collision bifurcation of a fixed point of a piecewise- 
smooth, continuous map are described by a piecewise-affme approximation of the form (J2]). 
The border-collision bifurcation occurs when the parameter, /i, is zero. As is varied from 
zero, invariant sets such as periodic solutions and topological circles may emanate from the 
bifurcation. The focus of this paper has been to fix fi at some nonzero value, assume ([2]) varies 
with other independent parameters, and investigate the nature of regions in parameter space 
where periodic solutions exist. In two-dimensional parameter space these resonance tongues 
commonly display a lens-chain structure. This paper has performed a rigorous unfolding 
about points for which resonance tongues have zero width - shrinking points. 

Curves corresponding to four different border-collision fold bifurcations form resonance 
tongue boundaries near non-terminating shrinking points. If {xi} denotes the corresponding 
5-cycle when it is well-defined, each boundary corresponds to the collision of one of the 
points x , x_d, ^(j-i)d and xu with the switching manifold. By Lem. @J at each boundary 
the corresponding border-collision matrix, (one of Ps, P$(-d), P s ai-i)d) and Pg(id)) is singular. 

At a shrinking point there exists an invariant polygon, V. If N = 2, or the shrinking point 
is terminating, V is planar, but in general V is nonplanar. V is comprised of uncountably 
many periodic solutions that have a given rotational symbol sequence, S. Under variation 
of parameters, V may persist as an invariant topological circle. In this case there exists 
curves along which one intersection point of the circle with the switching manifold maps to 
the other intersection point in a fixed, finite number of map iterations. We hypothesize that 
shrinking points generically occur densely along such curves, see [T5] . 

The two coexisting periodic solutions in a lens shaped resonance region have associated 
rotational symbol sequences with identical values for m and n but with values of I that differ 
by one. By adding or subtracting one from the value of / of these symbol sequences, one is 
able to obtain the symbol sequences of the periodic solutions in an adjoining lens. In short, 
throughout a lens-chain, the rotational number m/n is constant and / differs by one between 
lenses. Numerically we have observed that as n — > oo, the overall width of the lens-chain 
tends to zero and the number of individual lenses increases. We speculate that in the limit, 
a resonance tongue for a quasiperiodic solution of fixed frequency is a curve along which the 
fraction of points of the solution that lie in each half plane changes continuously. 

Shrinking points require solutions to intersect the switching manifold at two distinct 
points. For this reason shrinking points do not occur for ([2]) in one dimension because here 
the switching manifold is a single point. The one-dimensional, piecewise-linear, circle map 
studied in [15] is able to exhibit lens-chain structures because it has two switching manifolds. 

Resonance tongues corresponding to periodic solutions with non-rotational symbol se- 
quences need not exhibit a lens chain structure [IB]- But we believe that the collection of all 
symbol sequences that may generically exhibit a lens-chain structure is some class of which 
rotational symbol sequences are only the simplest type. For instance a period-n cycle of 
([2]) may undergo a Neimark-Sacker-like bifurcation when a complex conjugate pair of asso- 
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ciated multipliers crosses the unit circle. If this crossing occurs at A = e~^~ , where p, q G Z 
are coprime, then a period-ng orbit may arise in a manner akin to a terminating shrinking 
point. Numerically we have observed that such an orbit may exhibit a lens-chain shaped 
resonance tongue, even though its corresponding symbol sequence is non-rotational. This 
scenario corresponds to an n th iterate map of the same form as ([2D exhibiting lens-chains. 

The analysis presented in this paper applies to periodic solutions near the border-collision 
bifurcation of a fixed point. In general a periodic solution may result from global dynamics 
and enter regions of phase space bounded by many distinct switching manifolds. A study of 
the co dimension-two points resulting from such a periodic solution simultaneously colliding 
with different switching manifolds is left for future investigations. 



A Proofs of Lems. m and \5 



Proof of Lem. [7} We first show that (/ — A\_) and (J — A R ) cannot both be singular. Suppose 
otherwise. Then g T (I - A L ) = det(J - A L )ej = and g T (I - A R ) = det(J - A R )eJ = 0. 
Thus g T A L = g T A R = g T . Therefore g J M s = g T , i.e. M s has an eigenvalue, A = 1, which 
contradicts the assumption that (/ — Ms) is nonsingular. Thus at least one of (/ — A\_) and 
(/ — A R ) is nonsingular. Suppose w.l.o.g. that (/ — A\_) is nonsingular. 

For any b G M. N , let k = let c = b — kb and let y = fi(I — A\_)~ l c. Notice g T c = 0, 
thus the first component of y is zero, hence y = fic + Aiy = fic + A R y. Thus y = fiPsc + M$y 
and therefore y = — Ms)~ 1 Psc. 

Let Xq = kx + y. Multiplication on the left by ej yields s = ks . Using above results 
we have Xq = fi(I — Mg)~ 1 Pgb. Thus if b G nul^P^), then x = 0. Alternatively if b is chosen 
such that g T b = 0, then b = c and therefore xq = y. Hence if b G null(Ps) and g T b = 0, then 
y = and therefore c = b = 0. 

Now suppose Ps is singular. Then there exists b G null(P,s) with 6^0 and by the 
previous result, g T b ^ 0. Thus we have s = and k ^ 0, therefore s = 0, i.e. x lies on the 
switching manifold. 

Conversely suppose s = 0. Then s = — M s )~ 1 P s b = for any b G M N . Hence 

Ps is singular. □ 

Proof of Lem. {5[ Clearly if (J — M$) is nonsingular, (Tl3|) has the unique solution (fl6|) . To 
prove the converse, suppose now that (J — M$) is singular. We will show that (Tl3|) has no 
solution. Let K = Span(^) ± (the collection of all vectors orthogonal to g). Choose any 
c G K. We will show that the linear system, (/ — A L )x = fie, always has a solution, x, 
with s = ejx = 0. This is clear if (/ — A\_) is nonsingular for then x = — v4l) -1 c and 
s = faMzrAA = 0- Instead, suppose (/ — A\_) is singular. We now characterize the nullspace 
of (/ - A L j. 

We first show that if vq G null(J — A\_) and vo ^ 0, then ejv ^ 0. Suppose for a 
contradiction ejvo = 0. Let By denote the (N — 1) x (N — 1) matrix formed by removing 
the i th row and j th column from (J — A\_). Let v G IR^ -1 denote the last (N — 1) elements of 
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v . Then v 7^ and B^v = for each i, thus det(Bn) = for each i. In other words the first 
column of the cofactor matrix of (/ — A]_) is zero, equivalently g T = 0, which contradicts the 
assumption: g T b ^ 0. Hence ejv ^ 0. 

Let u, v G null(J — Then the linear combination, uejv — veju, is an element of the 
null space of (J — A\_) and the first element of this vector is zero. By the previous argument 
this vector must be the zero vector, thus u and v are linearly dependent. Hence null(J — A\_) 
is one- dimensional, therefore null(J — A\_) = Span(t> ). 

We have g J (I - A L ) = det(J - A L )eJ = 0, thus null((J - A L ) T ) = Span(g) = K ± . By 
the fundamental theorem of linear algebra, range (/ — A\_) = null((J — A]_) T ) ± = K. Notice 
fic G K, therefore there exists w G M. N such that (I — A\_)w = fic. Let x = w — ^-vo- 
Then (J — A\_)x = fic and s = 0. Summarizing, regardless of whether (J — A\_) is singular or 
nonsingular, we have found a vector, x, with zero first component, such that (I — A\_)x = fic. 

The vector x satisfies x = \ic + A L x = fic + A R x, therefore x = /iP^c + M s x and 
hence [iPsc G range(J — Ms). But c is arbitrary, hence Ps-fT C range(J — Ms)- K is 
(N — l)-dimensional thus P$K = range (J — Ms). By assumption, b K, thus f|T3l) has no 
solution. □ 



B Proof of Lem. ffl 

Proof. We begin with the case of non-terminating shrinking points. 

a) First, P§ = is singular by assumption, thus by Lem. 01 t — 0. Second, P^id) = 
PsWd) = f S ((i-i)d)o = Psca-ijd) is singular by assumption, where we have utilized Lem.[6bl 
for the second equality and (JTTj) for the last equality. Thus, by Lem. HJ tu = 0. 

b) Suppose for a contradiction, td = 0. Since also tid = from (a), by a double application 
of Lem. [H p solves the n-cycle solution system of S dld . From the definition, ffl9|) . it is 
seen that S dld = S^ d \ Therefore, pd solves the n-cycle solution system of S, hence 
p = pd- The point p is admissible by assumption, therefore p = pkd for any integer k. 
Putting k = m we obtain p = pi, hence p is a fixed point of and lies on the switching 
manifold. However cannot have a fixed point on the switching manifold because 
by assumption /j ^ and g T b ^ 0, so we have a contradiction. Thus td ^ 0. The 
remaining three points may be proven nonzero in a similar fashion. The given signs 
follow immediately from admissibility. 

c) Suppose for a contradiction, {p^ is of period n < n. Clearly h divides n. Then pm+k = 
Pk for any integers i and k. But from part (b), pd = pm+d lies left of the switching 
manifold. Thus Sm + d = L for i = 0, — 1, since the orbit is admissible. But 
gcd(d, h) = 1, hence in view of f|T^]) . we must have I > 2 + (" — l)n = n — h + 2 > f + 2. 
Similarly, also from part (b), = pm-d lies right of the switching manifold and hence 
Sin-d — R for i = 0, . . . , " — 1. Here it follows I < h < | which provides a contradiction. 
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d) From part (a), t = 0, thus by Lem. [TJ p solves the n-cycle solution system of S° = S. 
Similarly, since = 0, p solves the n-cycle solution system of S ld = S^~ d ' by fTiTJl) and 
therefore pd solves the n-cycle solution system of S. Let w{r) = rp + (1 — r)p d . Since 
{pi} is assumed to be admissible, {u>j(r)} is an admissible iS-cycle whenever < r < 1. 

From part (c), the n points {p{\ are distinct. We now show that the union of all 
the {wi(r)}, call it V, has no self-intersections. Suppose for a contradiction that 
w h( T i) — w i 2 ( T 2) f° r some %\ ^ z 2 and < ti,t 2 < 1. Let x = Wq(ti), then {x^ is an 
admissible 5-cycle and an admissible <S^ 2_ll )-cycle. Since r% ^ 0, 1, by part (b), so < 
and S-d > 0. Therefore S^-^ = L and S^-^-d = R. By Lem. 00, i 2 — %\ = 0, which 
is a contradiction. Hence V has no self-intersections and is therefore an invariant, 
nonplanar n-gon. 

To relate ([2]) to a map on the unit circle we describe a bijection between S 1 and V. The 
angular coordinate, 9 G [0,27r), uniquely describes a point on S 1 . Let j{6) = and 
T (6) = g -j(6). Let z : S 1 ^ P be defined by z(0) = w j{e)d (r(9)). It is easily verified 
the function z is a bijection. The induced map on S 1 is g(8) = (z^ 1 ofoz)(9) = # + 
i.e. rigid rotation with rotation number m/n. 



□ 



Proof. We now consider the case of terminating shrinking points. 

Let v = y + iz be an eigenvector corresponding to the multiplier A 



matrix A^. Let D 



' 2irm N 



n ' 
' 2irm ' 



sin 



cos 



2irm "> 



n ' 
' 2nm N 



Then Tfi 



cos( 
— sin 



2n jm \ 



e n for the 

2irjm \ 



n I 
2irjm \ 



sm, 

( 2-irjm \ 
COS( —z — ) 



and 



cos 
— sin 

A J L [y z] = [y z\T)i . Notice also D n = I. The first component of v must be nonzero because 
otherwise A R [y z] = [y z]D and hence Mg[y z] = [y z]D n = [y z], violating the assumption 
that (J — Mg) is nonsingular. Thus E c intersects the switching manifold. Furthermore, E c 
is two-dimensional, because otherwise there exists another eigenvector for A = e ± ^ R , call 



u, vejv 



vejv, is an eigenvector 



it v, linearly independent to v, and the linear combination, 
with zero first component contradicting the previous argument. 

Let x = xq = ay+ (3z + x*( L > for some as yet undetermined scalars a and (3. Let X{ denote 
the i th iterate of x via the symbol sequence <S = L n . Then Xi = [y z\D l [a (3} T + x*^ and hence 
{x^ is a period-n cycle. We wish to choose a and (3 such that so = ej[y z}[a P] T + s*^ = 
and s_d = ej [y z]D~ d [a P] T + 3*^ = 0. Combining these two equations yields the linear sys- 



tem, X[a p} T = -s*^ [1 1] T where X 



e{y 



ejz 



has determinant, det(X) 
It follows that 



_ e]ycos(f ) + ejzsm(^ ) -ejy sin(f ) + eJ"*cos(f ) 



((ejy) 2 + (ejz) 2 ) sin(^) ^ 0. Thus we let [a (3} T = -s^X^l 1] T . 



Sid 



*(L) 



1 



COS 



27r(i+i). 



COS 



(56) 



Thus the period-n cycle, is admissible and solves the n-cycle solution system of S. 

Therefore x = p. Thus {p{\ has period n which verifies part (c) of the lemma. Also ti — s», 
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for each i, thus by (|56|) we have parts (a) and (b). To verify part (d), note that as in the 
non-terminating case, p and pd are admissible solutions to the n-cycle solution system of S. 
Thus w(t) = rp + (1 — r)p d is an admissible solution to the n-cycle solution system of S 
for all r G [0,1). The union of all such cycles is an invariant, V, which lies on E c , so is 
planar. The points, Pi, lie on an ellipse and ordered so that V has no self- intersections and 
is an n-gon. As before, a bijection between S 1 and V may be constructed to show that the 
restriction of (T3J) to V is homeomorphic to rigid rotation with rotation number m/n. □ 
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